Mon. Not. R. Astron. Soc. 0000,[THT3l(2009) Printed 15 December 2009 (MN MbX style file v2.2) 



OS 

o 
o 

(N 
O 

Q 

in 



Azimuthally Symmetric Theory 0/ Gravitation (I) 

On the Perihelion Precession of Planetary Orbits 

G. G. Nyambuya* 

Received: 5 Sept. 2009 / Accepted with Moderate Revision: 9 Dec. 2009 



o 

>V 
Oh. 



> 

On 
(N 

<N 

On 
O 



- 1—1 

x 

S3 



ABSTRACT 

From a purely none-general relativistic standpoint, we solve the empty space Poisson equation 
(V 2 $ = 0) for an azimuthally symmetric setting, i.e., for a spinning gravitational system like 
the Sun. We seek the general solution of the form $ = $(r, 9). This general solution is con- 
strained such that in the zeroth order approximation it reduces to Newton's well known inverse 
square law of gravitation. For this general solution, it is seen that it has implications on the 
orbits of test bodies in the gravitational field of this spinning body. We show that to second 
order approximation, this azimuthally symmetric gravitational field is capable of explaining at 
least two things (1) the observed perihelion shift of solar planets (2) that the mean Earth-Sun 
distance mus t be increasing - this resonates with the observ ations of two independent groups of 
astronomers dKrasinskv & Brumberg Il2004t IStand ish 2005) who have meas ured that the me an 
Earth-Sun distan ce must be increasing at a rate of about 7.0 ± 0.2 m/cy dStandishll2005l) to 
15.0 ± 0.3 m/cy dKrasinskv & Brumberg 1120041) . In-principle, we are able to explain this result 
as a consequence of loss of orbital angular momentum - this loss of orbital angular momentum 
is a direct prediction of the theory. Further, we show that the theory is able t o explain at a satis - 
factory level the observed secular increase Earth Year (1.70 ± 0.05 ms/vr; iMiura et at 1 120091) . 
Furthermore, we show that the theory makes a significant and testable prediction to the effect 
that the period of the solar spin must be decreasing at a rate of at least 8.00 ± 2.00 s/cy. 

Key words: astronomical unit, azimuthal symmetry, orbit, perihelion shift, solar spin 



1 INTRODUCTION 

From as way back as the 1850s, it has been known that the or- 
bit of the planet Mercury exhibits a peculiar motion of its perihe- 
lion, specifically, the perihelion of Mercury advances by 43.1 ± 0.5 
arcsec/century. When Newton's theory of gravitation is applied to 
try and explain this (by making use of the oblateness of the plan- 
ets because when the Sun's gravitational force acts on the oblate- 
planets, the oblateness causes torque [on the planets] and this torque 
is thought to give rise to the anomalous mo tion of the plan ets); it 
was found first by Leverrier in 1859 see e.g. iKenvonl dl990} ) that it 
predicted a p recession of 53 2 arcsec/century which is larger than 
the observed jKenvonlll990l) . With the failure of Newton's theory to 
explain this, it was proposed that a small undetected planet was the 
cause. Careful scrutiny of the terrestrial heavens by telescopes and 
spaces probes reveals no such object - the meaning of which is that 
the cause may very well be a hitherto unknown gravitational phe- 
nomena - Einstein was to demonstrate that this was the case, that 
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there existed a hitherto unknown gravitational phenomena that is the 
cause of this peculiar motion. 

With the herald of Einstein's General Theory of Relativity (GTR) in 
1915, Einstein immediately applied his GTR to this problem; much 
to his elation which caused him heart palpitations - he obtained the 
unprecedented value of 43.0 arcsec/century and this was (and is 
still) hailed as one of the greatest triumphs for the GTR and this lead 
to its quick acceptance. Venus, the Earth, and other planets show 
such peculiar motion of their perihelion. Observations reveal a shift 
of 8.40 ± 4.80 and 5.00 ± 1.00 arcsec/century respectively (see 
e.g. Kenyon 1994). Einstein's theory is able to explain the perihelion 
shift of the other planets well, so much that it is now a well accepted 
paradigm that the perihelion shift of planetary orbits is a general 
relativistic phenomena. 

Einstein's GTR explains the perihelion shift of planetary orbits as a 
result of the curvature of spacetime around the Sun. It does not take 
into account the spin of the Sun and at the same time it assumes all 
the planets lay on the same plane. The assumption that the planets 
lay on the same plane is in the GTR solution only taken as a first or- 
der approximation - in reality, planets do not lay on the same plane. 
In this reading we set forth what we believe is a new paradigm; we 
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have coined this paradigm the Azimuthally Symmetric Theory of 
Gravitation (ASTG) and this is derived from Poisson's well accepted 
equation for empty space - namely V 2 $ = 0. Poisson's Law is a 
differential form of Newton's Law of Gravitation. We explain the 
perihelion shift of the orbits of planets as a consequence of the spin 
of the Sun - i.e. solar spin. It is well known that the Sun does exhibit 
some spin angular momentum - specifically, it [the Sun] undergoes 
differential rotation. On the averag e, it spins on its sp in axis about 
once in every ~ 25.38 days (see e.g. lMiura et a/.l20fM) . Its spin axis 
makes an angle of about 83° with the ecliptic plane. It is important 
that we state clearly here that by no means have we discovered a new 
theory or a set of new equations; we have merely applied Poisson's 
well known azimuthally symmetric solution to gravity for a spinning 
gravitating body. 

Further, with regard to Einstein's GTR - vis; in its solution to the 
problem of the perihelion shift of planetary orbits, it [the GTR] 
assumes the traditional Newtonian gravitational potential, namely: 
$(r) = ~GM/r, where G = 6.667 x lO -11 ^ - ^ -2 is New- 
ton's universal constant of gravitation, A4 is the mass of the central 
gravitating body and r is the radial distance from this gravitating 
body. Einstein's GTR which is embodied in Einstein's law of gravi- 
tation, namely: 




Figure 1. This figure shows a generic spherical coordinate system, with the 
radial coordinate denoted by r, the zenith (the angle from the North Pole; the 
colatitude) denoted by 6, and the azimuth (the angle in the equatorial plane; 
the longitude) by ip. 



(1) 



is designed such that in the low energy limit and low spacetime 
curvature such as in the Solar System, this equation reduces di- 
rectly to Poison's equation. In Einstein's law above, i? M „ is the 
Ricci tensor, 7? the Ricci scalar, g M „ the metric of spacetime, A is 
Einstein's controversial cosmological constant which at best can be 
taken to be zero unless one is making computations of a cosmolog- 
ical nature where darkenergy is involved, and k = 87rG/c 4 where 
c = 2.99792458 x 10 8 ms _1 is the speed of light in vacuum; and 
Poisson's equation is given by: 



V 2 $ 



4-irGp, 



(2) 



where p is the density of matter and the operator V 2 written for 
spherical coordinate system (see figure[T]for the coordinate setup) is 
given by: 



the confines of its error margins " - what immediately comes to mind 
is that can a theory have error margins or is it not experiments that 
have error margins? As will be seen, certain undetermined constants 
(A^) in the theory emerge and at present, one has to infer these from 
observations and it is here that the error margins of the ASTG come 
into play. 

Further, we show, that in-principle, the ASTG does explain (1) the 
increase in the mean Earth-Sun distance, (2) the increase in the mean 
Earth-Moon distance etc, and these emerge as a consequence of the 
fact that from the ASTG, the orbital angular momentum is not a 
conserved quantity as is the case in Newtonian's gravitational theory 
and Einstein's GTR. That the orbital angular momentum is not a 
conserved quantity may lead one to think that the ASTG violets the 
Law of Conservation of angular momentum - no, this is not the case. 
The lost angular momentum is transferred to the spin of the orbiting 
body and as well as the Sun. 
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As already been said, our solution or paradigm, hails directly from 
Poisson's equation, which in itself is a first order approximate so- 
lution to Einstein's GTR, albeit with the important difference that 
we have taken into account solar spin. This fact that our paradigm 
explains reasonably well - within the confines of its error margins; 
the precession of planetary orbits as a consequence of solar spin and 
at the sametime the GTR explains this same phenomena well as a 
consequence of the curvature of spacetime raises the question "Is 
the precession of the perihelion of solar orbits a result of (1) so- 
lar spin or (2) is it a result of the curvature of spacetime?" If 
anything, this is the question that this reading seems to raise. An an- 
swer to it, will only come once the meaning of the ASTG is fully 
understood. 

In the above we say the ASTG "explains reasonably well - within 



2 THEORY 

For empty space: V 2 "l? = 0; and for a spherically symmetric setting 
we have <!> = $(r) and this leads directly to Newtonian gravitation. 
For a scenario or setting that exhibits azimuthal symmetry such as a 
spinning gravitating body as the Sun we must have: $ = "l>(r, 0), we 
thus shall solve the Poisson equation: V 2 $(r, 9) = 0. The Poisson 
equation for this setting is readily soluble and its solution can read- 
ily be found in most of the good textbooks of electrodynamics and 
quantum mechanics for example - it is instructive that we present 
this solution here. 

We shall solve Poisson's equation for empty space (V 2 $ = 0) 
exactly; by means of separation of variables, i.e. we shall set: 
$(r, 9) = 4>(r)$(0). Inserting this into the Poisson equation we 
will have after some basic algebraic operations: 
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$(r) dr 



, d$(r) 
dr 



+ 



8 



$(6l) sinS 96> 



(4) 



The radial and the angular portions of this equation must equal some 
constant since they are independent of each other. Following tradi- 
tion, we must set: 



$(r) dr 



,d<f>(r) 



= £(£+1), 



and the solution to this is: 



$i(r) = Air + 



Br 



(5) 



(6) 



where At, and Br are constants and I — 0, 1, 2, 3, ... . If we set the 
boundary conditions; <&e(r = oo) = 0, then Ai — for all I. Now, 
just as Einstein demanded of his GTR to reduce to the well known 
Poisson equation in the low energy regime of minute curvature, we 
must demand that <f>(r), in its zeroth order approximation - where 
£ — and the terms I 1 are so small that they can be neglected; 
the theory must reduce to Newton's inverse square law; for this to be 
so, we must have: 



Bp = —Xtc 



(7) 




Figure 2. The elliptical planetary orbits have the Sun at one focus. As the 
planets describe their orbits, their major axes slowly rotate about the Sun 
in the process shifting the line from the Sun to the perihelion through an 
angle Aip during each orbit. This shift is referred to as the precession of the 
perihelion. 



where Xe is an infinite set of dimensionless parameters such that 
Ao = 1 and the rest of the parameters \i for £ > 1 will take values 
different from unity and these constants will have - for now, until 
such a time that we are able to deduce them directly from theory; to 
be determined from the experience of observations. In the discussion 
section, we shall hint at our current thinking on the nature of these 
constants. With this given, it means we will have: 



Mr) 



-XeC 



(GMV 



Now, moving onto the angular part, we will have: 



(8) 



(9) 



and a solution to this is a little complicated; it is given by the spher- 
ical harmonic function: 



$(6») = P £ (cos( 



(10) 



of degree I and Pi (cos 8) is associated Legendre polynomial. As al- 
ready said, the derivation of &(r, 8) just presented can be found in 
most good standard textbooks of quantum mechanics and classical 
electrodynamics. Since equation (J9} is a second order differential 
equation, one would naturally expect there to exist two independent 
solutions for every £. It so happens that the other solutions give infin- 
ity at 8 = (0, 7r), which is physically meaningless (see e.g. Grifitts 
2008). Now, putting all the things together, the most general solution 



is given: 



*(r,0) = -£ 



. 2 (GM\ e+1 p( f 



(11) 



which is a linear combination of all the solutions for £. In the case of 
ordinary bodies such as the Sun, the higher orders terms [i.e. £ > 1: 
of the term (GA4 /rc 2 ) ], will be small and in these cases, the 
gravitational field will tend to Newton's gravitational theory. Equa- 
tion dlll l is the embodiment of the ASTG, and from this, we shall 
show that one is able to explain the precession of the perihelion of 
planetary orbits. 

In this equation [i.e., il 1H. nowhere does the value of the Sun's spin 
(Tq ~ 25.38 days) enter into our equation. This may lead one to 
asking "So where has this been taken into account?". To answer 
this, it is important to note that if the potential is a function of r 
only i.e., $ = $(r), then, it technically is a function of r and 8 as 
well (with the ^-dependence being trivial). What this means is that 
spherical symmetry implies an azimuthal symmetry around any arbi- 
trarily chosen axis. If a specific axis is singled out, e.g. , by the spin of 
a body about the spin axis, then, the spherical symmetry of the static 
body is broken, and only an azimuthal symmetry remains and this 
azimuthal symmetry is only about the plane cutting the body into 
hemispheres such that this plane is normal to the spin axis. For any 
other plane cutting the body into hemispheres, the two hemispheres 
are asymmetric. From this we see that the azimuthally symmetric 
solution is consequence of the breaking of the spherical symmetry 
by the introduction of a spin axis, hence thus one is automatically 
lead to consider the solutions for which $ = <J?(r, 8). In this way, 
the spin has been taken into account. 



2. 1 Equations of Motion 

We shall derive here the equations of motion for the azimuthally 
symmetric gravitational field, $(r, 8). We know that the force per 
unit mass [or the acceleration i.e., g = — V"J>(r, 8)] is given by 



© 2009 RAS, MNRAS 0000,[TlfT3l 



4 G. G. Nyambuya 



a = (r — rip )f + (ri/5 + 2rtp)0 (see any good textbook on Clas- 
sical Mechanics) where a single dot represents the time derivative 
d/dt and likewise a double dot presents the second time derivative 
d 2 /dt 2 . Comparison of a = (f — rtp 2 )i + (rip + 2rip)6 with (g); 
i.e.: a = g, leads to the equations: 



dV 
dt 2 



dp 

df~ 



d$ 

dr ' 



(12) 



for the f-component and for the 0-component we will have: 
dr dip _ 1 d$ 



d 2 tp 



dt 2 



dt dt 



(13) 



d 2 u 
dtp 2 

and: 

dJ _ 
~dt ' 



1 d${u,6 p ] 
V d(9 D 



d$(u,g p ) 
d6> H ' 



dti 



1 d$(u,d p ) 



du 



(18) 



(19) 



This ends our derivation of the equations of motion for the field 
$(r, 9). Before we proceed to our main task of showing how equa- 
tions dl8land|T9l> explain the precession of planetary orbits, let us - 
for instructive purposes, first lay down Einstein's solution. 



Now, taking equation U3t and dividing throughout by rip and re- 
membering that the specific angular momentum J — r 2 tp, we will 
have: 



1 dip 2 dr 
tp dt r dt 



hence thus: 



l_d$ 

J dip 



ldJ 

l~dt 



l_d$ 

" J d0 ' 



dJ 
dt 



d$ 
" ~d9 ' 



(14) 



(15) 



The specific orbital angular momentum is the orbital angular mo- 
mentum per unit mass and unless otherwise specified, we shall refer 
to it as angular momentum. 

Digressing a little: what the above equation d 1 5 b means is that the 
orbital angular momentum of a planet around the Sun is not a con- 
served quantity. If it is not conserved, then the sum of the orbital 
and spin angular momentum must be a conserved quantity (if this 
angular momentum is not say transfered to the Sun or other solar 
bodies), the meaning of which is that at the different r-positions, the 
spin of a planet about its own axis must vary. This could mean the 
length of the day must vary depending on the radial position away 
from the Sun. We shall come to this later, all we simple want to do 
is to underline this, as it points to the possibility of a secular change 
in the mean length of the day. 

Now moving on, if we make the transformation u = 1/r, then for f 
and r we will have: 



dr T du , d r 
— — —J-r- and — — 
dt dip dt 2 - 



dJ du 

dt dip 



T 2 2 

J u 



dip 2 



(16) 



respectively. Inserting these into d!2t and then dividing the resul- 
tant equation by —it 2 J and remembering d!5t and also that dr = 
—du/u , one is lead to: 



d 2 u 
d^ 2 



1 d$(u,6)\ du 



■J- 1 



dip J dip 



+ u 



1 d<f>(u,8) 
I 2 du ' 



(17) 



The solutions that we shall consider are those for which 9 is a time 
constant, i.e. r — r(ip) and for the convenience we shall write 6 with 
subscript p, i.e., 9 P . This is just to remind us that it (9) is not a vari- 
able in the equations of motion as this is a constant for a particular 
planet p, hence: 



3 EINSTEIN'S SOLUTION 

When Einstein applied his newly discovered GTR to the problem 
of the precession of the perihelion of the planet mercury he obtained 
that the trajectory of solar planets must be described by the equation: 



d 2 u 
dtp'- 



+ u 



GM 3GMu 2 



J 2 



(20) 



where again u — 1/r. To obtain a solution to this equation, we note 
that the left hand side is the usual Newtonian equation for the orbit 
of planets, i.e.: 



d 2 u 
dp- 



+ u 



GM 
J 2 



(21) 



and the solution to this equation is: u = (1 + e cos tp) jl where e is 
the eccentricity of the orbit and I = (1 — e 2 )1Z where 1Z is half the 
size of the major axis of the ellipse. Written in different form, this 
solution is: 



1 



1 + e cos ( 



(22) 



where lZ m in is the planet's distance of closest approach to the Sun 
[see figure l[2j for an illustration]. This solution is a good approxi- 
mate solution to J20b because the orbit of Mercury is nearly Newto- 
nian. Consequently, we can rewrite the small term on the right hand 
side of (HI as: 3GM(l + ecos </j) 2 /Z 2 c 2 ; and in so doing, we make 
an entirely negligible error. With this substitution J20b becomes: 



d 2 u 



+ u 



GM 3GM 



dtp 2 J 2 ' l 2 c 2 

and the solution to this equation is 



1 + 2e cos tp + e cos tp\ 



1 + e cos tp 3GM 
u = ; h 



/ 



l 2 c 2 



1 + Y + 



6 



+ ttp sin p 



(23) 



,(24) 



Of the additional terms, the first i.e. (1 + e 2 /2) is a constant and the 
second oscillates through two cycles on each orbit; both these terms 
are immeasurably small. However, the last term increases steadily 
in amplitude with tp, and hence with time, whilst oscillating through 
one cycle per orbit; clearly this term is responsible for the precession 
of the perihelion. Dropping all unimportant terms we will have: 
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1 + e cos tp + erup sin ip 
I ' 



(25) 



where r\ — 3GM/lc 2 is extremely small. Thus all this leads us to: 

1 + ecos (Pew) 



I 



(26) 



where: (3b = (1 — v)- At the perihelion we will have: /3e<P = 2nir 
and this implies: ip = 2niif3g 1 ~ 2nir + 6nivGM/lc 2 . Essentially, 
this means that the perihelion advances by Aip — GtvGM /lc 2 per 
revolution and the resultant equation for the orbit is: 



1 + e cos (ip + Atp) ' 
hence thus the rate of precession of the perihelion is given by: 

Q-kGM 



\ T I E rc 2 (l- e 2 )TZ 



(27) 



(28) 



This is Einstein's formula derived in 1916 soon after he discovered 
the GTR. He [Einstein] concluded in the reading containing this for- 
mula: 

"Calculation gives for the planet Mercury a rotation of the orbit of 43" per 
century, corresponding exactly to the astronomical observation (Leverrier); 
for the astronomers have discovered in the motion of the perihelion of this 
planet, after allowing for disturbances by the other planets, an inexplicable 
remainder of this magnitude. " 



(3GM\ (GMV(3 



J J 



! cos 2 9„ 



(34) 



The left hand side of this equation (i.e. I32t is what one gets from 
pure Newtonian theory and the term on the right is the new term 
due to the first order term in the corrected Newtonian potential and 
likewise the second term on the right is a new term due to the second 
order term in the corrected Newtonian potential. 

Now, taking the term (3\u in equation ( 132b to the right hand side, we 
will have: 



d 2 u GM 2 

— + {l- Pl )u- — =p 2 lu . 



(35) 



We know that the solution of the right hand side of the above equa- 
tion when set to zero, i.e.: 



d 2 u 



GM 



is given by: 



/ 



1 + e cos(r]itp) ' 
where: 



(36) 



(37) 



fGMY ^ 2A lC os6> p 



(38) 



4 SOLUTION FROM THE ASTG 

For the present, we shall take the second order approximation of the 
potential $(r, 9) in-order to make our calculation for the precession 
of the perihelion of planetary orbits and this potential has been writ- 
ten down in l |29t . As has already been said; we shall consider only 
those solutions for which 8 is a time constant, i.e. r = r{tp) and 
for the convenience that we do not think of 9 as a variable we have 
set 9 := 9 P . The solutions r — r{tp) are those solutions for which 
the orbit of a planet stays put in the same #-plane. Now from the 
potential ( |29l we shall have: 



To obtain a solution to D5I > to first order approximation, we note that 
the left hand side has solution J37l l and that for nearly Newtonian 
orbits this solution u = (1 + ecos<p)/l, is a good approximation 
to d35 b for nearly Newtonian orbits such as Mercury for example. 
Consequently, we can rewrite the small term on the right hand side 



of d35 b as: 3G M ( 1 + e cos ip ) /I ; and make an entirely negligible 
error (see e.g. lKenvonlll990l) . With this substitution, equation J35 b 
becomes: 



dtp- 



2 

+ Vi u ■ 



GM 
~lf 2 ~ 



I 



1 



2e cos ip + e 2 cos 2 ip>) 



(39) 



dJ_ 
~dt 



^ GMu 



Ai sin Op + A2 



'3GM sin 261. 



2rc 2 



(31) 



Now making the transformation r = 1 /it, the first term on the left 
hand side of equation d30b transforms to: 



d 2 u GM 2 

dip z J 2 



where: 



and: 



(32) 



(33) 



and the solution to this equation is: 

1 + e cos rjiip (3-2 



- + ■ 



1 + 



e cos 2ip 



+ eip sin ip 



•(40) 



As before, i.e., as in the steps leading to Einstein's solution; of the 
additional terms, the first is a constant and the second oscillates 
through two cycles on each orbit; both these terms are immeasur- 
ably small. However, the last term increases steadily in amplitude 
with ip, and hence with time, whilst oscillating through one cycle 
per orbit; clearly this term is responsible for the precession of the 
perihelion. Now, dropping all the unimportant terms one is lead to: 



l + ecos77i<^ + e?72V?sin77ii^ 



(41) 
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<5>(u,9) = -GMu 



1 + Ai 



J \ du 



d 2 u 



dip 2 \ J 2 u 2 I dip 



( GMu 
\ c 2 

-GMu 2 



. . (GMu 
cos + A 2 I j — 

/ 2GMu cos 9 



1 + Ai 



3cos^ 



- A 2 



/3G.Mu\ 2 /3cos 2 0- 1 
I - 1? 



(29) 



,(30) 



where for the convenience we have set 7/2 = P2 and this quantity is 
extremely small, in which case cos 772 <p — 1 and sin 772 <p — 772 <p 
and using these approximations (in the cosine addition formula 

cos r\\ip + 772V sin 771 tp ~ cos 772 tp cos 771^3 + sin 772V? sin r]i<p = 
cos [(771 + 772)^]), we will have: 



1 + ecos [(771 + 772 )<£>] 



(42) 



Now, at the perihelion we are going to have: (771 + 7727V = 2n7r 
where n = 1, 2, 3, . . . and this implies tp — 2itn (771 + r/2) -1 = 
27rn[Vl-/3i + ~ 27rn[l - (/3i - 2 ( 3 2 )/2]- 1 = 27m[l + 

(/3i/2 — /3 2 ) + ...] hence: 95 ~ 2-7T71 + nX\h\ + nXihi, where: 



hi = 



and: 



7rG7V( \ / cos 9. 



2 ") 



(43) 



3 cos 2 6 V - 1 \ / 6?rGX \ 



12tt 



V ^ 2 / 



(44) 



This shows that per every revolution, the perihelion advances by: 



Aip 

T 



K\ ( 



lc 2 



3 cos 9 V 



12tt 



6tyGM \ 2 



lc 2 



and this can be written more neatly and conveniently as: 



Aip 



A<p 



cos 6 P I Aip \ 3 cos 2 9 P — 1 
— — Ai - ( ) — — — — — A 2 



\ T 1 1 



12-KT- 1 



•(45) 



This formula - which is a second order approximation; tells us of 
the perihelion shift of the planets. In the next section we will use 
this to deduce an estimate of the values of Ai and A 2 for the Solar 
System and thereafter proceed to calculate the predicted values of 
the perihelion shift. As a way of showing that these are solar values, 
let us denote (Ai and A 2 ) as (A® and A®) respectively. 



5 AN ESTIMATE FOR A® AND X§ 

If & p is the precession per century of the perihelion of planet 75, i.e. : 



then equation ( 145b can be written as: 



» p = £?p\f + gSpXf , 



where: 



(46) 



(47) 



and: 



T 1 1 



cos 9 r 



3 cos 



127T7— 1 



(48) 



(49) 



Given a set of the observed values for the size (l p ), the period of 
revolution r v , the tilt (9 V ) and the known precessional values of the 
perihelion of planets (£Pp bs ); these values are listed in columns 2, 3, 
4 and 8 of table (U( respectively; we can solve for A® and A° since 
2P V , stfp and SS V will all be known, thus one simple has to solve 
equation d47| > for any pair of planets as a simultaneous equation. 

The values of srf v and 3S P for all the solar planets are listed in 
columns 6 and 7 of table (Ul respectively. It is important that we 
state that the values of the Inclination listed in column 4 of table ^ 
are the inclination of the planetary orbits relative to the ecliptic plane 
and in-order to compute the inclination of these orbits relative to the 
solar equator we have to add 7° to this because the ecliptic plane 
and the solar equator are subtended at this angle. The solar equator 
is here defined as the plane cutting the Sun into hemispheres and this 
plane is normal to the spin axis of the Sun. 

Now, having calculated the values of Xf and A®, we will have to use 
these values (Xf and A®) to check what are the predictions for the 
precession of the perihelion of the other seven planets. If the predic- 
tions of our theory are in agreement with the observed precession of 
the perihelion of these seven planets, then our theory is correct and 
if the predictions are otherwise then, our theory cannot be correct - 
it must be wrong! 



For the present, we have calculated Xf and Af for the different 
planet pairs were we have all the information to do so and these 
values are displayed in table QO. The final adopted values are: 

Xf = 24.0 ± 7.0 and Xf = -0.200 ± 0.100. (50) 

and these values are the mean and the standard deviation - there is 
a 27% error in Xf and about twice (50%) that error margin in Xf . 
From the values given in d50t . the predicted values of the precession 
of the perihelion of the other seven planets i.e. Earth, Mars, Pluto; 
where computed and are listed in column 10 of table <UJ . The equiv- 
alent predictions of these values from Einstein's theory are listed in 
column 9 of the same table. Inspection of the predictions of our the- 
ory reveals that our predicted values are - as Einstein's predictions; 
in good agreement with observations. We believe that this does not 
mean the theory is correct but merely that it contains an element of 
truth in it. It means we have a reason to believe in it and as well a 
reason to peruse it further from the present exploration to its furthest 
reaches if this were at all possible! 

The reader should take note that in our derivation, we have assumed 
as a first order approximation the Newtonian result namely that the 



© 2009 RAS. MNRAS 0000, mfHI 



On the Perihelion Precession of Planetary Orbits 

Perihelion Precession of Solar Planetary Orbits According to the ASGT 



7 



Planet 



(AU) (yrs) 



Wind. 

(°) 



Precession (l"/100yrs) 



n 



Mercury 


0.39 


0.24 


7.0 





206 


3.50 


X 


10" 




1.72 


X 


10 2 




43.1000 ± 0.5000< c > 


13 


50000 


12 


80000 ±0.10000 


Venus 


0.72 


0.62 


3.4 





007 


5.19 


X 


io- 


1 


2.88 


X 


10 1 




8.0000 ± 5.0000< c > 


8 


62000 


12 


00000 ± 3.00000 


Earth 


1.00 


1.00 


0.0 





017 


1.57 


X 


io- 


1 


3.80 


X 


io- 


1 


5.0000 ± 1.0000^ 


3 


87000 


4 


00000 ± 1.00000 


Mars 


1.52 


1.88 


1.9 





093 


7.02 


X 


io- 


2 


2.43 


X 


io- 


2 


1.3624 ± 0.0005^ 


1 


36000 


1 


70000 ± 0.50000 


Jupiter 


5.20 


11.86 


1.3 





048 


3.02 


X 


10" 


3 


1.00 


X 


io- 


5 


0.0700 ± 0.0040( £ ) 





06280 





07000 ± 0.02000 


Saturn 


9.54 


29.46 


2.5 





056 


7.59 


X 


10" 


4 


1.72 


X 


io- 


7 


0.0140 ±0.0020^ 





01380 





01900 ± 0.00050 


Uranus 


19.2 


84.10 


0.8 





046 


1.09 


X 


io- 


4 


9.76 


X 


io- 


5 


(/) 





00240 





00250 ± 0.00070 


Neptune 


30.1 


164.80 


1.8 





009 


3.98 


X 


io- 


5 


9.13 


X 


io- 


11 


(/) 





00078 





00270 ± 0.00070 


Pluto^ 


39.4 


247.70 


17.2 





250 


5.77 


X 


10- 


5 


9.48 


X 


io- 


12 


(/) 





00042 





00140 ± 0.00040 



Notes : At the 2006 annual meeting of the International Astronomical Union, it was democratically decided that the solar 
test body Pluto is not a planet but a dwarf planet. For our purpose, its inclusion here as a planet is not affected by this 
decision for as long as this test body orbits the Sun like other planets. ^ The values of l v ,r p , Inc. and Ecc. are adapted from 
ISaganl ll974l) . ( c ) Adapted fr om lKenvonl ll990l) . (*» Adapted from |PitievaN2005l) . (e) Obtained by adding the extra precession 
determined by |Pitievdj200ll) and found in llorioU2008bl) to the standard Einsteinian perihelion precession. W Because of their 
long orbital duration covering at least 2 human lifetimes, no data is currently available covering one full orbital revolution for 
Neptune and Pluto hence there is not yet any observational values for the precession of their perihelia. The data for Uranus is 
unreliable (see e.g. Iloriol2008bl) . 

Table I. Above, column 1 gives the name of the planet p, column 2 gives l v which is the observed value for the orbital size 
of planet p, column 3 gives t p which is the period of revolution of the planet for one full orbit, column 4 is the tilt 9 P in 
degrees of the planet's orbit orbit relative to the ecliptic plane, column 5 gives the eccentricity of the orbit of the planet, while 
columns 6 and column 7 give the computed values jz/ p and SS V and column 8,9 and 10 give (1) the observed, (2) the GTR 
and (3) the ASTG precessional values of the planet. 



Estimation of the Values xf and 



Planet Plair 


Ai 


A 2 


Mercury-Venus 


15.8 


-0.0716 


Mercury-Earth 


32.8 


-0.4174 


Mercury-Mars 


20.0 


-0.1574 


Mercury- Jupiter 


26.1 


-0.2895 


Mercury-Saturn 


27.6 


-0.3112 


Mean 


:24.0 


-0.200 


Standard Deviation 


: 7.0 


-0.100 


Percentage Error 


:27% 


50% 



Table II. Column 1 in the first part of the table gives the name of the pair 
of the planets which have been used to obtain the pair of A-values listed 
in columns 2 and 3. In the second part of the table we compute the Mean, 
Standard Deviation and Percentage Error of the A-values 



angular momentum is a time constant. From the preceding section, 
clearly this is not the case. We have only assumed this as starting 
point of our exploration. It is hoped that taking into account the fact 
arising from the ASTG that orbital angular momentum is not a con- 
served quantity should lead to improved results that hopefully come 
closer to the observed values. 



6 NONE CONSERVED ORBITAL ANGULAR 
MOMENTUM AND ITS IMPLICATIONS 

Through equation J31b which clearly states that the orbital angular 
momentum of a planet must change with time; three immediate con- 
sequences of this are (1) a change in the mean Sun-Planet distance 
(2) a changing length of a planet's day and (3) a secular change in 
solar spin. In the subsequent subsection, we shall go through these 
implied phenomena. 



6.1 Increase in Mean Sun-Planet Distance 

One of the most accurately determined physical parameters in as- 
tronomy is the mean Earth-Sun distance which is about the size of 
the As tronomical Un it (AU) where 1 AU = 149597870696.1 ± 
0.1m i Pitieval 120051) and this is known to an accuracy of 10 cm 
dPitieval l2005h . The Astronomical Unit according to the Inter- 
national Astronomical Union (Resolution No. 10 197fQ) is de- 
fined as the radius of an unperturbed circular orbit that a mass- 
less body would revolve about the Sun in 27r/fc days where k — 
01720209895AU 3/2 day' 1 is Gauss' constant. This definition is 
such that there is an equivalence between the AU and the mass of 



1 see http://www.iau.org/static/resolutions/IAU1976 French.pdf 
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the Sun M which is given by GMq = k 2 A A . So, if Mq is fixed, 
it is technically incorrect to speak of a change AU. 

Before it was noticed that the mean Earth-Sun distance was chang- 
ing it made perfect sense to refer to the mean Earth-Sun distance as 
the Astronomical Unit. Now, (1) because units must not change, and 
(2) because of this fact that the mean Earth-Sun distance is chang- 
ing; then, until such a time that the Astronomical Unit is correctly 
defined so that it is a true constant as physical unit must be, it makes 
sense only to talk of the mean Earth-Sun distance instead of the As- 
tronomical Unit. 

That the m ean Earth-Sun distance is cha ngin g, this has bee n 



measured by Krasinskv & Brumberg ] J2004I) and Istandishl J2005h . 
iKrasinskv & Brumberg I J2004I) finds 15.0 ± 4.0 m/cy wh ich in 
SI units is (4.75 ± 1.27) x 10" 9 m/s and Istandishl feOOSh finds 
7.00 ± 0.20 m/cy which in SI units is (2.22 ± 0.06) x 10" 9 m/s 
where Icy — 100 yr. 

To this rather surprising res ult, i.e . , the ap parent secular change in 
the mean Earth-Sun distance. lloriol J2008ah states that the secular in- 
crease in the mean Earth-Sun distance can not be explained within 
the realm of classical physics. Contrary to this, we believe and hold 
that the ASTG can in-principle explain this result. The ASTG is 
well within the provinces of classical physics hence thus this re- 
sult is explainable from withi n the domain s and confines of clas- 
sical physics. In his reading dlorioll2008d) argues that the Dvali- 
Gabadadze-Porrati braneworld scenario - a none-classical theory, 
which is a multi-dimensional model of gravity aimed to the ex- 
planation of the observed cosmic acceleration without darkenergy, 
predicts, among other things, a perihelion secular shift, due to Lue- 
Starkman Effect of 5 x W 4 arcsec/cy for all the planets of the 
Solar System. It yields a variation of about 6m /cy for the increase 
in mean Earth-Sun distance; this is compatible with the observed 
time rate of change of the mean Earth-Sun distance hence giving the 
Dvali-Gabadadze-Porrati braneworld theory some breath. 

Ilorid d2008ah goes on to say that the recently measured corrections 
to the secular motions of the perihelia of the inner planets of the 
Solar System are in agreement with the predicted value of the Lue- 
Starkman effect for Mercury, Mars and, at a slightly worse level, 
the Earth. We shall show that in-principle, the ASTG can explain 
this result as a consequence of the none-conservation of orbital an- 
gular momentum of planets in this azimuthally symmetric gravita- 
tional setting. The none-conversation of the orbital angular momen- 
tum leads directly to a time variation in the eccentricity of planetary 
orbits. This makes the secular change a purely classical result. 

Now, given the definition of the eccentricity of an orbit: 



(TZmin 
TZmax 



(51) 



where IZmin and lZ ma x are the spatial extent of the minor and major 
axis respectively; and then, differentiating this with respect to time, 
one is lead to: 



de 
'dt 



K 

nun nun 

R 

in in 

tin 

max 



7Z 2 

' v in 



dt 



dt 



(52) 



dt 



and from this it follows that: 



de _ ( Tim. 

€ di ~ ~ { w, 



(7 + 1) 



iz„ 



(dTZ in i n 



(53) 



(54) 



and multiplying by H m in both sides, and thereafter substituting 
TLminj'TLma.x on the right hand side we will have: 



elZ nn n 



de 



dt 

therefore: 

dTZ in i n. 



-(i-e 2 )(i-( 7 + i) vTT^)) 



dTZ 

mi n 

dt 



dt 



(l-e3)(l-( 7 + l)VT=p) 



, (55) 



(56) 



Now, on the average, the time change of the minor axis must to a 
large extend be a good measure of the time change of the average 
distance IJZ) between the planet and the Sun, hence thus: 



dt 



(l-e 2 ){( 1+ l)VT 



de\ 



(57) 



In the realm of Newtonian gravitation where spherical symmetry is 
assumed thus producing equations only dependent on the radial dis- 
tance r, the eccentricity is an absolute time constant, i.e. de/dt = 0, 
and this directly leads to d IJZ) /dt = 0, hence when one finds that 
the mean Earth-Sun distance is increasing, it comes more as a sur- 
prise. If we consider azimuthally symmetry in Poisson's equation as 
has been done here, the result emerges naturally because the eccen- 
tricity is expected to increase with the passage of time - this we shall 
demonstrate very soon. 

In §(4), against the clear message from the ASTG, we assumed that 
the orbital angular momentum of a planet is a conserved quantity. It 
turns out that taking this into account leads us to two type of orbits 
(1) spiral orbits (2) the normal elliptical orbits with the important 
difference that the eccentricity of these orbits varies with time and it 
is this variation of eccentricity that we believe the secular increase 
of the mean Earth-Sun distance is rooted. 

Doing the right thing and taking into account the predicted change 
in the angular momentum, then equation ( |35t will be: 



d 2 u 
dp 2 



1 dj\ du 



J 2 u 2 dt J dip 



2 



GM a 

—p~ = P2lU , 



(58) 



and taking the change of angular momentum to first order approxi- 
mation from equation fl31b . one will have: 



dJ_ 

It 



. (GMY . 
Ai I I sm9 f 



-Actu 



(59) 



where a is clearly defined from this equation i.e. : 



There is no reason to assume that the rate of change of the minor and 
major axis be the same, thus we must set: 



°=2 



sin % 



(60) 
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it therefore follows that: 



d u 2a du 2 GM 



dtp* ltd^ +nlU ~ — =P2lU 



J 2 



and writing k — a/J 2 which is: 



(GM\ 



Aj [GM\ . 



(61) 



(62) 



where the Newtonian approximation J 2 = GMl has been used and 
K = GM I J 2 , the above becomes: 



d u du 2 2 
- k- — h 77! it — K = /J 2 lu . 
dip 1 dip 



(63) 



If the orbital angular momentum varies constantly with time, then 
J — Jt + Jo where Jo is the angular momentum at time t = and J 
is a time constant, then k — k(t) and K — K(t), meaning k(t) and 
K(t) will dependent not on the coordinates r, 9, ip but only on time, 
hence in solving the above equation we can treat these as constants 
since they do not dependent on r, 9, p. We believe the assumption 
that J = constant is justified because if that was not the case, there 
could be an accelerated increase in the orbital angular momentum 
and this could have been noticed by now. In this assumption that 
J = constant, we must have J being so small that it is not easily 
noticeable as it appears to be the case since we have had to relay 
on delicate observations to deduce the secular increase of the mean 
Earth-Sun distance. 

Now, to obtain a solution to this equation (i.e. 163b . we need first to 
get a solution to: 



d 2 u du 2 GM _ 
__ 2fc _ +r7lW __ = , 



and to obtain a solution to this, we need first to solve: 



d u du o n 



(64) 



(65) 



and to its solution we add GM / J 2 . The axillary differential equa- 
tion to this differential equation is: X 2 — 2kX + n\ = and its (i.e. 
equation[65} solutions are: 



X = k± \k 2 ~n\ = k±in 3 , 



(66) 



where 773 = \Jn\ — k 2 . If (773) 2 < the solution is: u = 
A e (k+va)v> _j_ 5 e (fc-t)3)*> w here A and B are constants, thus adding 
GM/J 2 we have: 



(67) 



u = Ae {k+m)v + Be {k -^ )v + 

J- 

and if (773) 2 = the solution is: u = (Aip + B)e kv thus adding 
GM/J 2 we have: 



u= {Aip + B)e KV + 



^ , GM 
J 2 ' 



(68) 



The solutions {67} and ( 168b are clearly spiral orbits. These solutions 
are obvious very interesting but because our focus is not on them, 
but on the solutions giving elliptical orbits in which the eccentricity 
varies, we shall not be looking into these spiral orbit solutions any 
further than we have already done. 

Now, in the event that (?73) 2 > the solution to equation l !64b is: 



1 + ee v cos(na<p) 
I ' 



(69) 



Now, using the same strategy as that used in §(3) and (4) to solving 
equations J20b and l !35b respectively, one finds that the resultant orbit 
equation will be: 



/ 



1 + ee fc ^ cos[(t72 + rj:j)ip] ' 



(70) 



and as before, at the perihelion we will have (772 +773)^ = 2-™ and 
this implies tp = 27177(772 + ?73) _1 — 27T7i[/32 + \/vt ~ k 2 ]" 1 = 
27rn[/3 2 + s/l - ft - fc 2 ]" 1 ~ 2tto[1 + (2/3 2 -/3i)/2- k 2 /2}' 1 
and taking only first order terms we will have: ip ~ 27rn[l + (/3i — 
2/3 2 )/2 + k 2 /2] and this shows that the perihelion will precess by 
an amount A(p = 27r[(/3i/2 — /3 2 ) + k 2 /2], and in comparison 
with Aip ~ 2n[f3i/2 — /3 2 ] obtained without taking into account 
the changing angular momentum, there is an additional precession 
of (Aip)+ ~ irk 2 . The value of k 2 for the Solar System is so small 
that in practice, one can neglect it, thus, we have not missed out 
much in our calculation in which we have assumed a constant orbital 
angular momentum. While this result is important our main thrust is 
to deduce the variation of the eccentricity of elliptical orbits (we 
shall shelf any deliberations on this result for a further reading). 

In equation ( 170b . the term te kip in the denominator is the eccentric- 
ity, let us write this as e* = ee k(p , and from this we see that the 
eccentricity varies with time - i.e.; as the orbital angular momen- 
tum changes with the passage of time, so does the eccentricity. Now 
plucking this into equation f57l > we can determine the variation of the 
mean Earth-Sun distance if we have knowledge of 7, unfortunately 
we do not have this. However, if we are to reproduce the observed 
variation of the Earth-Sun distance, one finds that if they were to set 
773 = 1.48 x 10~ 4 , which practically means that the orbit grows 
evenly at every point, one is able to explain the secular increase of 
the mean Earth-Sun distance. 

It should be said that, if the ASTG is to stand on its own - i.e., 
independent of observations, then it must be able to explain the result 
7b = 1.48 x 10~ 4 from within its own provinces. It is for this reason 
that we say, in-principle, the ASTG is able to explain the secular 
increase in the mean Earth-Sun distance and only until such a time 
when one is able to derive say the value 773 = 1.48 x 10 -4 from 
within the theory itself, will we be able to say the ASTG explains 
the secular increase in the Earth-Sun distance. 

Other than the secular increase in the mean Earth-Sun distance, there 
is also the in crease in the mean Earth-M oon distance. This has been 
measured bv lWilliams & Bo ggs (2009) to be ~ 3.50 x 10~ 3 m/cy 
and in SI units this is 1.11 x 10 -12 m/s. This observation provides 
a test for the ASGT, but unfortunately, we do not have the value of 
Ai so as to check what the ASTG says about this. We believe one 
cannot use the same A's values obtained for the Sun because these 
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values must be specific to the gravitating body and may very well 
be connected to the spin or the gravitating body in question. We are 
working on these ideas to improve the ASTG and at present we can 
only say it is prudent to assume that the A-values are specific to the 
body in question hence one has to calculate them from observational 
data. For the Earth, this increase in the Earth-Moon distance is but 
the only observations we have in-order for us to deduce Xf hence the 
ASTG is unable to make any predictions on this as it stands in the 
present. We hope in the future one will be able to deduce a general 
form of the Af-values, thus placing the ASTG on a level where it is 
able to make predictions that are independent from observations. 

Important to note from e» = ee v is that, as ip i — > — oo, the ec- 
centricity will decrease and the reserve is that the eccentricity will 
increase as ip i — > +00 decreases. An increasing eccentricity leads 
to a secular decrease in the Planet-Sun distance and a decreasing ec- 
centricity leads to a secular increase in the Planet-Sun distance. This 
means the sense in wh ich the planet orbits the Sun is i mpo rtant! Be- 
cause we believe from lKrasinskv & Brumberg] d2004l) and lStandishl 
( l2005h . that there is a secular increase in the Earth-Sun distance, this 
means the current direction of rotation of the Earth around the Sun 
must be such that ip 1 — > —00. This must be true for other planets 
rotating in the same sense as the Earth; and to any (object in the So- 
lar System) that rotates in the direction opposite to this, this body 
will experience a secular decrease in its distance from the Sun. 



6.2 Secular Increase in the Orbital Period of Planets 

Given that through the passage of time - what is suppose to be a 
sacrosanct parameter - the mean Earth-Sun distance; is changing, 
and that the time change of the specific orbital angular momentum 
is given j = 2rr6+r 2 6, then, if as in the case of Newtonian gravita- 
tion the specific orbital angular momentum of a planet is a conserved 
quantity, i.e. J = 2rr6+r 2 6 = 0, then accompanying this result of a 
changing mean Earth-Sun distance must be an increase in the length 
of a planet's duration for one complete orbit since 6/6 = —2r/r. 
Given that 8 = 2n /Ty where Ty is the orbital period of a planet, 
the equation 6/6 = —2r/r becomes: Ty /Ty = 2r/r. Plucking in 
the relevant values for the Earth, one is lead to Ty — 2.97 ms/cy. 
Since T® — 365.25Td where T§ is the period of an Earth day, it 
follows that: = 365.257^, further, follows that we must have: 
T§ = 8.13 fis/cy - this value is at odds with physical reality; for 
records held for over 2700 yrs indicate that the Earth day changes 
by an amount "T® — +1.70 ± 0.05 ms/cy (see e.g. iMiura etal] 
120091) . which is about 200 times that expected if the orbital angular 
momentum where a conserved quantity as in Newtonian gravitation 
- clearly, this suggests that the orbital angular momentum may not 
be conserved. 

If say the conserved quantity where the total angular momentum of 
a planet, i.e. the sum total of the spin angular momentum (S) and 
the orbital angular momentum, then S = —J and if the radius of 
the planet is not changing with time, then Td = — 2-kTZ 2 JT^. For 
the Earth, one finds that T, 



D 



-5.18 s/cy which is ~ 3000 times 
the observed value - this can not be, sure something must be wrong. 
We shall explain this observational value T§ — 1.70±0.05 ms/cy 
from the ASTG. 

From the ASTG, we have: 



theory 



(6.00 ±2.00) x 10 _1 V 



and we know that: 



1) -Al 



7y_ 

Ty 



(71) 



(72) 



hence plucking in the observed values and remembering not to forget 



that for the Earth = 365. 257^, then we will have: 



= -(2.28 ± 0.07) x 10~ 15 s~\ 



(73) 



This value - vis, the order of magnitude, is on a satisfactory level in 
good agreement with observations. We take this as further indication 
that the ASTG contains in it, a grail of the truth. 



6.3 Secular Increase in Solar Spin 

We know that angular momentum must be conserved but according 
to J 3 1 b . it is not conserved. This lost orbital angular momentum must 
go somewhere - it cannot just disappear into the thin interstices of 
spacetime or into the wilderness of spacetime thereof. Let Ltot be 
the sum total angular momentum of the Solar System, were we con- 
sider that the Solar System is composed of the planets. If the sum 
total of the angular momentum of a planet and its system of satellite 
is Jp Dt , then C to t = MqSq + J2i M^J\ ot . We would expect that 
the total angular momentum of the Solar System be conversed, that 
is dCtot /dt = 0. From this we must have: 



Mq 

Mo S Q 



-y 



'Mi 




M e 





(74) 



and djp 0t /dt — dj p /dt hence thus: 



r 







27^0 M e 
T b TZ q M e 2-kTI 2 



o 



E 



Mi djj 
Me dt 



(75) 



and this means the orbital period of the Sun must be changing. If 
we assume that the Sun's radius has remained constant through the 
passage of time, i.e. 72-q = (which is certainly not true), then what 
we obtain from the above is a minimum value for the secular change 
in the Sun's spin. The reason for invoking this assumption is because 
there currently is no information o n the secular change of the Sun's 
radius (see e.g. IMiura ef"aZll2009l ). hence we make this assumption 
so that we can proceed with our calculation. As already said, what 
we get is not the exact secular change in the Sun's spin but a lower 
limit to this. 

The second term in equation {75}, i.e. Me /Me', represents the 
effect of solar mass loss, which can be evaluated in the follow- 
ing way. The Sun has a luminosity of at least 3.939 x 10 26 W, 
or 4.382 x 10 9 kg/ s; this inclu des electromagnetic radiation and 



the contribution from neutrinos l lNoerdlinge3l2008h . The particle 



mass loss rate due to the solar wind is ~ 1.374 x 10 kg/ s 
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(see e.g. iNoerdlin ger 2008). From this information, it follows that 
M G j/Mq ~ 9.10 x lO-^raT 1 . 

Now, the last term in equation J75t can be evaluated from the 
ASTG since J is known - so doing, one finds that it is equal to 
~ -(4.00±1.00)xl0 _6 cy _1 ; this implies T = 8. 00±2.00 s/cy. 
This result is a significant 10 6 times larger than the term emerg- 
ing from the solar mass loss so much that we can neglect this al- 
together and consider only the last term in equation J75b hence 
Tq = 8.00 ± 2.00 s/cy. This value i s sign ificantly larger com- 
pared to that calculated bv lMiura et al.\ d2009h where these authors 
find a value of 21.0ms/cy. Currently, no serious measurements on 
the secular change in the period of the solar spin has been made. It 
should be possible to u ndertake this effort and with respect to the 
ASTG, and the result of lMiura etal\ <2009h . this experiment would 
act an arbiter. 

Furthermore, the authors iMiura et al. propose that the Sun 

and the Earth are literally pushing each other away (leading to the 
increase in the Earth-Sun distance) due to their tidal interaction and 
they believe that this same process is what's gradually driving the 
moon's orbit outward: they say "Tides raised by the moon in our 
oceans are gradually transferring Earth's rotational energy to lu- 
nar motion. As a consequence, each year the moon's orbit expands 
by about 4 cm and Earth's rotation slows by about 30/is". Further 
IMiura efall d2009l) assumes that our planet's mass is raising a tiny 
but sustained tidal bulge in the Sun. They calculate that, thanks to 
Earth, the Sun's rotation rate is slowing by 30/j.s jcy. Thus according 
to their explanation, the distance between the Earth and Sun is grow- 
ing because the Sun is losing its angular momentum - the ASTG 
gives a different explanation altogether and this is in our opinion, 
very interesting. 



7 DISCUSSION AND CONCLUSIONS 

We have considered Poisson's equation for empty space and solved 
this for an azimuthally symmetric setting - we have coined the 
term Azimuthally Symmetric Theory of Gravitation (ASTG) for 
the emergent theory thereof. From the emergent solution, we have 
shown that the ASTG is capable of explaining certain observed (and 
yet to be observed) anomalies: 

(1) Precession of the perihelion of planets. 

(2) Secular increase in the Earth-Sun distance. 

(3) Secular increase in the Earth Year. 

(4) Secular decrease in solar spin. 

(5) Spiral orbits must exist. 

One of the draw-backs of the ASTG as it currently stands is that it is 
heavily dependent on observations; for the values of Xe need (have) 
to be determined from observations. Without knowledge of the X' e s, 
one is unable to produce the hard numbers required to make any 
quantifications. Clearly, a theory incapable of making any numerical 
quantifications is useless. This must be averted. We shall make use 
of the solar values of the \' e s in shading some light into our current 
thinking on this, i.e. finding a general form for the constants Xi; 
In the subsequent paragraphs, we shall make what we believe is a 
reasonable suggestion and give our current envisage-ment on the 
general form for these constants. 



(1) First things first, if the constants Xe where all independent of 
each other, then, the theory would clearly be horribly complicated. 
If we take as guide the Principle of Occam's Razor which in most if 
not all cases, leads to the simplest theory, then, these constants must 
be dependent on each other somehow so as to reduce the labyrinth 
of complications thereof. The simplest imaginable such dependence 
is Xi = F(£)Xi\ in this way, the entire system of constants Xe 
is dependent on just the one constant Ai. This idea that the sys- 
tem of constants be dependent on just one constant is drawn from 
the theory of polynomial functions where for a polynomial func- 
tion F(x) = "}2^L c n x n , one can have "well behaved" polynomial 
functions for which the constants c n have a general form, i.e., were 
they dependent on n; e.g., e x = IC^o 2 "™/ 71 - We eny i sa g e tne 
function $(r, 8) to be a "well behaved" function. By "well behaved" 
we simple mean its system of constants, Xe, is critically dependent 
on I just as the constants c n depend on n. 

(2) Second, we could like that on a practical level, only the second 
order approximation of the theory must suffice, this means the terms 
I > 3 must be practically negligible. We have already shown herein 
that the second order approximation of the ASTG is able to explain a 
sizable amount of anomalous observations. With the ASTG written 
in its second order approximation and as will be shown in the second 
reading (i.e., a follow-up reading that we hope will be published in 
the present journal), one is able without much difficulties, to explain 
from this second order approximation, the emergence of molecu- 
lar bipolar outflows in star forming systems, as a gravitational phe- 
nomena. If the other terms beyond the second order approximation 
become practically significant, one will have difficulties to explain 
outflows. So in a way, we are not going to pretend but clearly state 
that, we want - albeit with a priori and posteriori justification; to fine 
tune the theory so that it is able to explain the emergence of bipolar. 
This is the strongest reason we want the terms for which £ > 3 to be 
so small such that in practice one can neglect them entirely. 

(3) Third and most important, the only data point we have of these 
constants is the determined values for the Sun, i.e., Xf = 24.0 ± 7 
and A2 = —0.2 ± 0.1. If logic is to hold - as it must; then, our 
suggestion, Xe = F(i)Xi; must be able to explain this. We find that 
the following proposal: 

f(-l) e+1 \ 

Xe = [jfyuFjj Al ' (76) 

meets (1), (2) and (3). We shall assume this result until such a time 
evidence to the contrary is brought forth. Checking on (3) we see that 
within the error margins A© ~ [(-1) 2+1 / ((2 2 )!(2 2 )]A . Further 
checking on (2); from JTSJl we will have A 4 = 3.40 x 10~ 30 Ai 
which is practically small; the meaning of which is that all terms for 
which £ > 3 can in practice be neglected entirely. 

If the above proposal proves itself to be correct, then, the resultant 
theory will have just one undetermined parameter Ai . We are not go- 
ing to try and deduce what this parameter depends on but simple hint 
at our current thinking. We believe this parameter must depend on 
the angular frequency of the spin of the gravitating body in question. 
If we can find the correct dependence, then, the ASTG will stand on 
it own thus positioning itself on the podium to make testable predic- 
tions. We have left the task to make this deduction an exercise for 
the follow-up reading. 
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The fact the we have deduced the crucial parameters Af & Af from 
experience, means we have in the current reading done some reserve 
engineering. Normally, a theory must give these values and make 
clear predictions, just as when Einstein wrote down his equations 
and found that his theory predicted a factor 2 difference when com- 
pared to Newton's theory when it come to the bending of light by the 
Sun and when applied to the Sun-Mercury system, it accounted very 
well for the then unexplained 43.0" per century for the precession 
of the perihelion of the orbit of this planet; it just came out right. 
There were no free parameters that needed fitting as is the case of 
the ASTG. As argued above, once a general form for the Xi is found, 
this setback of the ASTG will be solved. Because we were able to 
obtain the values A & Af which lead acceptable values for the per- 
ihelion precession, means that the values srf & SB are not random but 
systematic. If the theory was all wrong, then, only luck would make 
the obtained values for si & 8$ give values of A® & A® such that 
equation J47l > give in general, acceptable values for the precession of 
the perihelion of the planets. 

With regard to the values obtained for the precession of the perihe- 
lion of solar planets, it can be said that, the values obtained from 
the ASTG as shown in column 10 (table U when weighed against 
the observational values listed in column 10 (of the same table) are 
acceptable. Given that we have taken into account the fact the or- 
bits of these planets are not found laying in the same plane, this can 
hardly be a coincident or an accident since changing their inclina- 
tion by just 1° will alter the predicted values of the precession of 
their perihelion. 

Ilorid d2008d) states that the secular increase in the mean Earth- 
Sun distance cannot be explained within the realm of classical 
physics. Contrary to this, we believe and hold that herein - we have 
shown from within the provinces of classical physics that this re- 
sult is explainable from within the domains and confines of classical 
physics. Before the present, the reason why perhaps this observation 
appeared beyond the reach of classical physics is because classi- 
cal physics has not really considered gravitation as an azimufhally 
symmetric phenomenon as has been done in present reading. This 
strongly suggests that the ideas presented herein need to be explored 
further for they contain a debris of the truth. 

One of the interesting outcome that was not explored in this reading 
for fear of digression is that the ASTG has a provision for spiral 
orbits (equation 1671 and [68). These orbits occur when (n^) 2 ^ 0. 
This condition implies the existence of a region (r $C lZ cr it) in 
which spiral orbits will occur. Evaluating the inequality (723) 2 $J 0, 
leads to: 1Z cr it = (2AiCM /c 2 ) cos 2 ((9/2), and from this, it is easy 
for one to deduce that spiral orbits are unlikely in the Solar System 
since these will have to occur inside the Sun because lZcrit 72-0 . 

At this point as we approach the end of this reading, we feel strongly 
that we must address the question; "Does the spin along the az- 
imuthal axis of a gravitating body induce an azimuthal symmetry 
into the gravitational field for this spinning body?" To answer this, 
we must ask the question; "Will a contracting none-spinning cloud 
of gas experience any bulge alone its equator?" First, we know that 
the equatorial bulge will occur on a plane perpendicular to the spin 
axis. Since a none-spinning gas cloud is going to have to spin axis, 
there is going to be no spin axis about which the equatorial bulge 
will occur. If the material in the cloud is randomly and uniformly 
distributed, the cloud will exhibit a spherically symmetric distribu- 
tion of mass and its gravitational field is expected to be spherically 



symmetric. A spherically spherically symmetric gravitational field 
is one that only has a radial dependence, i.e. ip = ip(r). 

Now, if the gas cloud is spinning, the centrifugal forces will cause 
there to exist a disk and the material distribution will have an az- 
imuthal symmetry, i.e. p = p(r, 9). Should not this azimuthal sym- 
metric distribution of matter induce an azimuthal gravitational field? 
From Poison's equation l[2}, p — p(r,9) implies <I> = $(r, 8); 
should not this, i.e. $ = 5>(r, 8), hold as-well for a body spinning 
gravitating body in a vacuum? From this, clearly, a spinning gravitat- 
ing body ought to exhibit an azimuthal symmetry. We have argued in 
the last paragraph of §(2) that the spin of a gravitating body breaks 
the existing spherical symmetry of the non-spinning body and the 
above argument is just adding more to this. It is from this that the 
subtitle and running head finds its justification. 

If the ASTG turns out to be correct - as we believe it will; then, we 
have an important question to ask; "What is the speed of light do- 
ing in a theory of gravitation because from ([7} we see that the con- 
stants G and c are intimately tied-up together ? This is a similar if 
not a c ongruent question that has been asked bv lMartin & Anderson! 
j2009t) in their expository work on Earth Flyby Anomalies (AFA). 
The empirical formula deduced to quantify EFA contains in it 
the s peed of light, c, so in their exposition of the phenomena of 
AFA. iMartin & AndersorJ fe009l) have asked the perdurable ques- 
tion "What is the speed of light doing there?". EFA are thought to be 
a gravitational phenomena, so, what does the speed of light have to 
do with gravitation - really? If there is an intimate relationship be- 
tween the speed of light and gravitation, then, one will be forgiven 
to think this suggests a link between gravitation and the theory of 
light - electromagnetism. The speed of light, c, appears to be dire to 
the ASTG presented herein. Why not another value but the speed of 
light, c? We shall leave these matters hanging in-limbo. 

In relation to the question above, i.e., "What is the speed of light do- 
ing in a theory of gravitation", one notes that Newtonian gravitation 
- which requires instantaneous interaction as a postulate; does not 
imply the dependence of the gravitational potential on the azimuthal 
angle for a spinning body as is the case in the ASTG, because at any 
instant t, the gravitating body appears spherically symmetric. Here 
we have the speed of light c coming in because of the azimuthal 
symmetry. Does this speed of light c link (or not) the propagation of 
the gravitational phenomena to the speed of light? In the present, we 
can only pause this as a question, for we still have to do further work 
on these ideas. 

In closing, allow me to say that we find it hard to call what has 
been presented herein as "A New Theory of Gravitation". When 
one tells you they have come up with a new theory of gravitation, 
what immediately comes to mind is that they have discovered a new 
principle upon which gravitation can further be understood from 
the present understanding. The ASTG is not founded on any new 
physical principle but on the well known vintage equation of Pois- 
son. What we have done is simply taken the azimufhally symmetric 
equations of this equation and applied them to gravitation. Based 
on this understanding, it is difficult to call it a new theory. Yes, the 
azimufhalfy symmetric equations of Poisson have brought new and 
exciting physics - perhaps only because of this, the title of this read- 
ing finds its qualification. 
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